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This paper analyses theoretically and numerically the effect of varying grating amplitude on the 
extremely asymmetrical scattering (EAS) of bulk and guided optical modes in non-uniform strip- 
like periodic Bragg arrays with stepwise and gradual variations in the grating amplitude across the 
array. A recently developed new approach based on allowance for the diffractional divergence of the 
scattered wave is used for this analysis. It is demonstrated that gradual variations in magnitude 
of the grating amplitude may change the pattern of EAS noticeably but not radically. On the 
other hand, phase variations in the grating may result in a radically new type of Bragg scattering — 
double-resonant EAS (DEAS). In this case, a combination of two strong simultaneous resonances 
(one with respect to frequency, and another with respect to the phase variation) is predicted to take 
place in non-uniform arrays with a step-like phase and gradual magnitude variations of the grating 
amplitude. The tolerances of EAS and DEAS to small gradual variations in the grating amplitude 
are determined. The main features of these types of scattering in non-uniform arrays are explained 
by the diffractional divergence of the scattered wave inside and outside the array. 



I. INTRODUCTION 



Extremely asymmetrical scattering (EAS) of waves in 
periodic arrays is a new type of Bragg scattering that is 
realized when the scattered wave propagates parallel or 
almost parallel to the boundaryfies) of a strip-like peri- 
odic array [llllllllEi,|iElalI3,ll3 This type 
of scattering is radically different from the conventional 
Bragg scattering in periodic arrays. For example, EAS is 
characterized by a strong resonant increase of the scat- 
tered wave amplitudes inside and outside the array; the 
smaller the grating amplitudCjlhe lar ger the amplitudes 
of the scattered waves H H |1 B S [lOi El- In addi- 
tion, the incident and scattered waves inside the array 
each split into three waves H 0, S H ElIJl • Two of 
these scattered waves and two of the incident waves in- 
side the array are evanescent waves which are localized 
near the array boundaries 0, Q, 111 • The third scattered 
wave is a plane wave propagating at a grazing angle into 
the array ^ i . 

Our recent publications 0, S 111 have demonstrated 
that the physical reason for EAS is related to the diffrac- 
tional divergence of the scattered wave inside and out- 
side the array. A new powerful approach for sim- 
ple analytical analysis of EAS, based on allowance for 
this diffractional divergence, was introduced and justi- 

fied iHiiiiaiia 

In the case of bulk electromag- 
netic waves, this approach was shown to give the same 
coupled mode equations as the dynamic theory of scatter- 
ing that waspreviously used for the theoretical analysis 
of EAS H HH S i, 0. However, the dynamic the- 



ory of scattering is not suitable for analysis of EAS of 
guided and surface waves on account of extremely awk- 
ward calculations involved, while the new approach is 
readily applicable for all types of waves, including sur- 
face and guided optical and acoustic modes ^J^^ J^] • 

EAS has enormous potential for new important practi- 
cal applications in the development of novel optical com- 
munication devices (e.g. narrow-band optical filters, res- 
onators, couplers, switches, lasers, etc), optical sensors 
and measurement techniques. For example, the strong 
resonant increase of the scattered wave amplitude during 
EAS may result in high-quality EAS-based resonators 
and high sensitivity of sensors and measurement tech- 
niques. New non-coUinear geometry of EAS can lead to 
the development of highly tuneable optical and ultrasonic 
devices, and may also result in an improved side-lobe 
structure of filtered signals. The possibility of concen- 
tration of the scattered wave energy in narrowchannels 
can be used for amplification and lateral compression of 
waves, as well as for effective coupling of a planar waveg- 
uide and a fibre. 

However, manufacturing EAS-based structures and de- 
vices will inevitably be related with various imperfections 
of periodic arrays. Because of the resonant character 
of EAS, it is possible to expect that at least some of 
these imperfections may be crucial for the experimen- 
tal observation and practical use of this type of scat- 
tering. Thus, a very important practical problem is to 
investigate theoretically the effect of array imperfections 
on EAS, and determine the tolerance of the scattering. 
In addition, non-uniform arrays are often manufactured 
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on purpose. For example, non-uniform chirped gratings 
can effectively compress pulses broadened due to dis- 
persion of optical modes in a slab or optical fibre 
(i.e. the dispersion can be compensated). Launching of 
solitons into non-uniform nonlinear Bragg gratings can 
be significantly easier than for uniform nonlinear Bragg 
gratings \13L ll^ . The side- lobe structure of filtered sig- 
nals can be noticeably improved (suppressed) by using 
non-uniform gratings with slowly varying grating ampli- 
tude 

There are also two other special reasons for using 
nonuniform arrays in the case of EAS. Firstly, edge effects 
at the array boundaries may result in noticeable undesir- 
able energy losses in the scattered wave, caused by an ad- 
ditional re-scattering of the scattered wave. These losses 
may be noticeable because of the large amplitude of the 
scattered wave propagating along the boundaries. Non- 
uniform arrays with gradually increasing grating ampli- 
tude can be used for a substantial reduction of this effect. 
Secondly, nonuniform arrays with varying grating ampli- 
tude may result in a radically new type of EAS , as 
well as cause concentration of the wave energy within 
narrow channels inside the array In addition, non- 
uniform arrays may result in a significant reduction of 
the relaxation time to steady-state EAS. 

Therefore, EAS of bulk or guided opticalwaves has 
been analysed for non-uniform arra ys with step-like 
variations of the grating amplitude |lCt ITU . It was 
demonstrated that a step-like variation in magnitude 
of the grating amplitude may result in a noticeable re- 
arrangement of the intensity distribution of the scattered 
wave inside the array [Tll |. The sensitivity of EAS to 
small uniform and non-uniform (step-like) variations in 
the grating amplitude was also determined ini. Step- 
like phase variations in the grating inside the array were 
shown to have a much stronger effect on EAS than the 
magnitude variations In thin arrays, they resulted 
in a radically new type of EAS double- resonant extremely 
asymmetrical scattering (DEAS) that is characterized by 
a unique combination of two sharp simultaneous reso- 
nances with respect to frequency and the phase variation 
in the grating ,10j. 

However, non-uniform arrays with step-like variations 
of the grating amplitude will hardly reduce edge effects 
at the array boundaries. Moreover, the presence of ad- 
ditional interfaces at which the step-like variations take 
place must increase undesirable edge effects and the re- 
sultant energy losses. This is especially the case for 
DEAS, where the scattered wave amplitude is exception- 
ally large 1^. Step-like variations in the grating ampli- 
tude will also hardly improve the side-lobe structure of 
processed (filtered) signals. Finally, manufacturing peri- 
odic arrays usually results in imperfections characterized 
by slow (but not step-like) variations in the grating am- 
plitude (e.g. due to non-uniform etching, photolithogra- 
phy, etc). However, EAS and DEAS in non- uniform ar- 
rays with slow variations of the grating amplitude have 
not been analysed previously, though such analysis is cru- 



cial for successful development of new practical applica- 
tions of these types of scattering. 

Therefore, the aim of this paper is to analyse theoret- 
ically and numerically the steady-state EAS and DEAS 
of bulk and guided electromagnetic waves in non- uniform 
periodic arrays with slowly varying grating amplitude. 
Coupled wave equations describing EAS and DEAS in 
arrays with slowly varying grating amplitude are derived. 
The field structure in the scattered waves is determined 
inside and outside the array by means of numerical so- 
lution of the coupled wave equations. The tolerance of 
EAS and DEAS to small slow variations in the grating 
amplitude (array imperfections) is also investigated the- 
oretically. The results are compared with those obtained 
for EAS and DEAS in uniform and non-uniform arrays 
with step-like variations of the grating amplitude. 



II. COUPLED WAVE EQUATIONS 

The structure under investigation and the geometry of 
scattering are presented in figure 1. An incident wave 
is scattered in a strip-like periodic array of width L. 
The period of the grating is assumed to be constant in 
the whole array, while the grating amplitude can vary 
in phase and/or in magnitude across the array, i.e. in 
the x-direction (figure 1). We will consider only stepwise 
and/or slow (i.e. small at distances of the order of the 
grating period) variations in the grating amplitude. In 
this case, the grating amplitude can be assumed to be 
locally constant at each point of the array, except for the 
values of the x-coordinate at which the steplike variations 
take place (at x — Li in figure 1). The array is assumed 
to be uniform along the y-axis, i.e. in the directions par- 
allel to the array boundaries. 




x = X - L'^ X = L 



FIG. 1: Scheme for EAS in non-uniform strip-like periodic 
Bragg arrays with varying grating amplitude and the total 
array width equal to L. 

In sections 2-5 we analyse EAS and DEAS of bulk TE 
optical waves in periodic Bragg arrays represented by a 



3 



periodic variation of the dielectric permittivity: 

Cs = e + ei{x) exp(iq • r) + e^(x) exp(-iq ■ r) 

if < x < L, 

es = e if a; < 0, or X > L, (1) 

where the niean dielectric permittivity e is the same in 
all parts of the structure (inside and outside the array), 
the amplitude of the grating ei(x) is small: 



|ei(a;)|/e < 1, 



(2) 



and q is the reciprocal lattice vector (|q| ~ 2tt/A, where 
A is the grating period). There is no dissipation of elec- 
tromagnetic waves inside or outside the array, i.e. e is 
real and positive. The complex grating amplitude ei(a;) 
varies in magnitude and/or phase within the array. 

A plane TE electromagnetic wave (with the electric 
field parallel to the z-axis) is incident on the array at an 
angle Oq (measured counterclockwise from the x-axis — 
figure 1). We assume that the Bragg condition is satisfied 
precisely: 



ki - ko 



(3) 



where ko is the wavevector of the incident wave, ki is 
parallel to the array boundaries (figure 1), |ki| = |ko| — 
ko — cje^/^/c, uj is the angular frequency, and c is the 
speed of light in vacuum. 

If condition (2) is satisfied for all values of x, then 
the amplitudes of the incident and scattered waves vary 
slowly inside the array, i.e. their variations at distances 
of about one wavelength are small as compared with the 
values of these amplitudes. In this case the approxima- 
tion of slowly varying amplitudes is valid and only two 
waves — incident and scattered — need to be taken into ac- 
count inside and outside the array: 

E{x) = Eq(x) exp{ikoj;X + ikoy-y — iLut}+ 

Ei{x)e:>ip{ikoy - iujt}, (4) 

where Eq(x) and Ei{x) are the slowly varying amplitudes 
of the electric field in the incident and scattered waves, 
respectively, kox — ko cos^Oiand koy_= kg sin^o- 

It was shown previously [1 S H E IH that in the 
geometry of EAS, there are two opposing mechanisms de- 
termining the behaviour of the scattered wave amplitude. 
On the one hand, the scattered wave amplitude must in- 
crease along the direction of its propagation (i.e. along 
the y-axis) due to scattering of the incident wave inside 
the array. On the other hand, the scatteredwave ampli- 
tude must decrease along the y-axis due to a significant 
diffractional divergence of this wave 9i 10, ij. 

The steady-state EAS occurs when the increment in the 
scattered wave amplitude caused by the scattering is 
exactly compensated by the decrement caused by the 
diffractional divergence 0, H S i, ElIS- The new 
approach for the theoretical analysis of EAS is based 



on the separate determination of each of the two con- 
tributions to the scattered wave amplitude. The cou- 
pled wave equations describing the steady-state EAS 
are then derived from the comparison of these contri- 
butions S^J^jSini- 
Here, we apply this approach to the considered case 
of EAS and DEAS in non-uniform periodic arrays with 
slowly varying grating amplitude. During the first step, 
we determine the contribution due to scattering, disre- 
garding the diffractional divergence. In this case the am- 
plitude of the scattered wave must increase along the 
y-axis. If the amplitudes of the incident and scattered 
waves vary slowly inside the array, then at any point in 
the array, the incident wave amplitude and the grating 
amplitude can be considered to be locally constant. Thus 
the scatteringinduced increments in the amplitude of the 
scattered wave along the direction of its propagation at 
any point of the array are determined only by the local 
values of the grating amplitude and the amplitude of the 
incident wave, regardless of the type of Bragg scatter- 
ing. Therefore, we can find these increments by means of 
the conventional dynamic theory of scattering in uniform 
arrays. 

The coupledwave equations in the conventional d^ 
namic theory of scattering are well known 

M. 



iventional dy- 

miiiiiQ 



d£'o/da;o 
d£'i/da;o 



iViEi, 



(5) 
(6) 



where Fq and Fi are the coupling coefficients, and the 
XQ-axis is parallel to the reciprocal vector of the grating 
(see figure 1). Equation (6) gives the scattering-induced 
rate of changing the scattered wave amplitude along the 
xg-axis. The rate of changing amplitude of the scattered 
wave along the direction of its propagation can then be 
obtained from equation (6) by means of the simple substi- 
tution of dxo = — cos(pi/2 — 77-|-6'o)dy = — siii{ri — 9o)dy: 



V dy 



-iToEoix)sm{7j-9o). (7) 



scattering 



where 77 is the angle (measured counterclockwise) be- 
tween the positive XQ-direction and the wavevector of the 
incident wave (figure 1). 

During the second step, we disregard the scattering 
and consider only the diffractional divergence of the scat- 
tered wave (beam). In this case, the rate of decreasing 
amplitude of the scattered wave can be derived by substi- 
tuting the entire scattered field Esc = Ei exp(ifciy — iujt) 
into the Helmholtz equation V^Esc + k\Esc = 0: 



dx'^ 



o2 o 

El + ^Ei -f 2ifc— £;i 
ay^ ay 



0. 



(8) 



Here, we can neglect the second-order derivative of the 
slowly varying amplitude Ei with respect to y as com- 
pared with the first-order derivative. This gives us the 
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parabolic equation of diffraction 



dy 



divergence 



2fci 9x2 



(9) 



that determines the rate of decreasing scattered wave am- 
phtude along the direction of its propagation due to the 
diffractional divergence. 

As has already been mentioned, in the steady-state 
EAS, the contributions to the scattered wave amplitude, 
caused by the scattering and diffractional divergence, 
must exactly compensate each other. Therefore, the sum 
of rates (7) and (9) must give zero. This condition results 
in the following equation: 



dx^ 



where 



Ko = -2fciro sm{rj-0o). 



(10) 



(11) 



Note that in equations (8)-(ll) we have deliberately 
used ki instead ko, even though ki — ko for bulk electro- 
magnetic waves. This was done to make the presented 
derivation valid for all types of waves, including surface 
and guided modes, for which ki may be not equal to ko 
(e.g. for scattering of TE modes guided by a slab into 
TM modes of the same slab). 

Equation (10) is one of the coupled wave equations 
describing EAS (or DEAS) in a non-uniform periodic 
array. The second equation is obtained from equation 
(5). Recall that in the considered geometry of scatter- 
ing (see figure 1), the incident wave amplitude must de- 
pend only on the x-coordinate. Therefore, substituting 
dxo = dxi COST? = dxcosrj/ cos^o into equation (5) gives: 



dEojx) 
dx 



iKiEiix), 



where 



Ki — Fi cos?// cos 00- 



(12) 



(13) 



Equations (10) and (12) arc the complete set of coupled 
wave equations for EAS (or DEAS) of electromagnetic 
waves in a non-uniform periodic array with varying grat- 
ing amplitude. As has already been mentioned above, 
these coupled wave equations are valid for the description 
of EAS and DEAS of all types of waves, including sur- 
face and guided optical modes. The difference between, 
for example, bulk and guided optical waves is only in dif- 
ferent values of the coupling coefficients Fq and Fi. For 
bulk electromagnetic waves 8] 



Fo -F* -eluj^/[2c^ko cos 77] 



(14) 



while for guided modes, these coefficients are determined 
using one of the conventional dynamic theories of scat- 
tering for optical slab modes [13, HI; uM (^^e section 
6 for more detail). 



There is no coupling between incident and scattered 
waves outside the array, i.e. at x < and x > L. There- 
fore, the coupling coefficients Ko,i are equal to zero out- 
side the array, and the coupled wave equations (10) and 
(12) take the form: 



d^Ei{x)/dx^ ^ 0; dEoix)/dx = 0. 



(15) 



One can easily see that the coupled wave equations 
(10), (12) and (15) in the case of non-uniform arrays with 
slowly varying grating amplitude appear to have exactly 
the same form as for uniform arrays [1,0 SB- However, 
for non-uniform arrays, the coupling coefficients Ki and 
Ko in equations (10) and (12) are dependent on the x- 
coordinate because the grating amplitude ei depends on 
X. This results in substantial difficulties with analytical 
solution of these equations. Therefore, in sections 3 and 
4, we present results of numerical analysis of the coupled 
wave equations (10) and (12) for periodic arrays with 
several typical dependences of the grating amplitude on 
the x-coordinate. The numerical method and boundary 
conditions used are described in the appendix. 



III. EAS IN NON-UNIFORM ARRAYS 

In this section, we analyse numerically the effect of 
slowly varying grating amplitude on EAS of bulk TE elec- 
tromagnetic waves in periodic Bragg arrays. The results 
obtained for the non-uniform arrays are compared with 
those for EAS in uniform arrays 0|. We also investi- 
gate numerically the tolerance of EAS to small gradual 
variations (imperfections) in the grating amplitude. 



x = x = L 

(i) constant 




x=0 x=L 

(iii) sinusoidal 



FIG. 2: Three different profiles of the grating amplitude in- 
side a periodic array: (i) constant grating amplitude (uniform 
array), (ii) linearly varying grating amplitude (non- uniform 
array), and (iii) sinusoidally varying grating amplitude (non- 
uniform array). The mean grating amplitudes are the same. 

For example, consider three different types of depen- 
dences of the grating amplitude on the x-coordinate in- 
side the array — figures 2(i)-(iii). In figure 2(i), the grat- 
ing amplitude is constant throughout the array, i.e. we 
consider a uniform array of width L and the grating am- 
plitude eiQ. EAS in such arrays was investigated previ- 
ously for bulk 8] and guided 6, "^l optical waves. Mathe- 
matically, the profile of the grating amplitude in such an 
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FIG. 3: The dependences of the relative scattered wave amphtudes on distance from the front array boundary inside the arrays 
with eio = 5 X 10~^, e = 5, 9o = tt/4, the wavelength in vacuum A = 1 jim (the grating has a period of 0.58 pm and is inclined 
at the angle of n/8 to the front array boundary); (i) L = 12 pm, (ii) L = 20iim, (iii) L — 40 pm, (iv) L = 80 pm. Curves (a): 
uniform array with constant grating amplitude eio — figure 2(i). Curves (b): non-uniform array with the linear dependence of 
ei(x) — figure 2(ii). Curves (c): non-uniform array with the sinusoidal dependence of ei(a;) — figure 2(iii). 



array is given by the equations: 
ei{x) = 



if < a; < L, 
otherwise. 



(16) 



The second of the considered arrays is a non-uniform 
array with a linear dependence of magnitude of the grat- 
ing amplitude on the x-coordinate in such a way that the 
grating amplitude is zero at the front and rear boundaries 
of the array (figure 2(ii)): 

{ieiox/L ifO<a:<L/2, 
ieio{L~x)/L ifi/2<a;<L, (17) 
otherwise. 

The gradient of changing grating amplitude in (17) is 
chosen such that the mean grating amplitude in the non- 
uniform array is the same as in the uniform array de- 
scribed by equations (16). Thus the integral 



ei{x)dx 



(18) 



must be the same for both the dependences presented 
by equations (16) and (17). This condition is necessary, 
because EAS is strongly dependent on the grating ampli- 
tude HiLilSlliSj. If, for example, ei{x) increased gradually 
from zero to eio, the mean grating amplitude would be 
noticeably smaller than that of the uniform array, result- 
ing in a significant increase in the scattered wave ampli- 
tude. If we assume that the mean grating amplitudes 
are the same for the uniform and nonuniform arrays of 
the same widths (i.e. integral (18) is the same for both 
the arrays), then we will be able to investigate only the 
effects of non- uniformity of the grating on EAS. 

The third of the considered arrays is characterized by 
sinusoidal variations of the grating amplitude so that the 
grating amplitude is again zero at the front and rear array 
boundaries (figure 2 (iii)): 

J eio[l - cos(27ra::/L)] if < x < L, 
£1 a;) = < , 19) 

^ ' lO otherwise. ^ ' 

Here, the mean grating amplitude is again the same as for 
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arrays (16) and (17), i.e. integral (18) is the same for all 
three dependences (16), (17) and (19). The relationship 
between the maxima of the grating amplitudes in arrays 
(16), (17) and (19) can be written as 



max{ei(a;)}iincar = max{ei(a;)}sinusoidai 



= 2 max{ei(x)} 



stepwise 



= 2ei 



(20) 



The selection of non- uniform arrays in the form of (17) 
and (19) is not arbitrary. Analysis of arrays with zero 
grating amplitude at the front and rear boundaries (as for 
arrays in figures 2(ii) and 2(iii)) is practically important 
because such arrays may significantly improve the side- 
lobe structure of filtered signals and substantially 
reduce unwanted edge effects at the array boundaries. 

Figure 3 presents the dependences of the scattered 
wave amplitudes on the x-coordinate inside periodic ar- 
rays (16), (17) and (19). It can be seen that for all array 
widths, slowly varying grating amplitude does not in- 
troduce radical changes in the pattern of scattering. It 
can also be seen that the effect of slowly varying grating 
amplitude on amplitudes of the scattered waves is notice- 
ably smaller for narrow arrays. For example, for arrays of 
12 pm thickness (figure 3(i)), the difference between the 
scattered wave amplitude Ei in the non-uniform array 
with sinusoidal variation of the grating amplitude (19) 
and the scattered wave amplitude in the uniform array 
(16) is less than 0.4%. For linear dependence (17) this dif- 
ference is even smaller (figure 3(i)). At the same time, it 
can reach 10-15% in arrays of 40 pm width (figure 3(iii)), 
and over 25% in arrays of 80pm width (figure 3(iv)). 

This can be explained by the diffractional divergence 
of the scattered wave. In narrow arrays, the diffractional 
divergence effectively smooths out variations in the am- 
plitude of the scattered wave, caused by different values 
of the grating amplitude. For such arrays, the scattered 
wave amplitudes are determined (to a very good accu- 
racy) by an average value of the grating amplitude (fig- 
ures 3(i) and (ii)). However, the diffractional divergence 
is effective only within a limited distance. This is be- 
cause when the scattered wave spreads along the x-axis 
(due to the diffractional divergence) from one part of the 
array to another, it experiences rescattering in the grat- 
ing. As a result, the wave can diverge along the x-axis 
only until it is re-scattered by the grating, i.e. it can 
spread only within a finite distance (see also |0|). This 
distance must increase with decreasing mean grating am- 
plitude, because the efficiency of re-scattering decreases 
with decreasing grating amplitude. It should also in- 
crease with increasing gradient dei (x) / dx, because in this 
case variations of the amplitude of the scattered wave 
along the wavefront will increase, resulting in stronger 
diffractional divergence. Typically, the distance through 
which the scattered wave can spread inside the array be- 
fore being re-scattered is about several where 
7 = (i^o^i)'/'. 

In our examples (see equations (16), (17), (19) and fig- 
ures 2 and 3), the critical array width is about ~ 37^^"'^, 



where 7av is the average value of 7 in the considered ar- 
rays. Therefore, if L < 37^;^^, then the scattered wave 
amplitude is very well determined by the average value 
of the grating amplitude inside the array (figures 3(i) and 
3(ii)). If the array width L > 37a^^, then the scattered 
wave amplitudes inside a non-uniform array are signifi- 
cantly affected by local values of the grating amplitude, 
and the pattern of scattering differs from that for the 
uniform array with the average grating amplitude (fig- 
ures 3(iii) and 3(iv)). 

So far, we have analysed non-uniform arrays with sig- 
nificant variations of the grating amplitude. Another 
practically important problem is the effect on EAS of 
small grating imperfections, e.g. variations in the grat- 
ing amplitude, which are inevitable during manufacture 
of periodic arrays. Therefore, in this paper, we analyse 
tolerance of EAS to small variations of the grating am- 
plitude across the array (i.e. in the x-direction). 




x = X- L x = x = L 

(i) constant (ii) linear 



x-O x-L 

(iii) linear 




X - x-L 

(iv) sinusoidal 




X = X-L 

(v) constant 



FIG. 4: Five different profiles of the grating amplitude with 
small variations inside a periodic array: (i) constant grating 
amplitude eio (uniform array), (ii) linearly increasing grating 
amplitude (non-uniform array), (iii) linearly decreasing grat- 
ing amplitude (non-uniform array), (iv) sinusoidally varying 
grating amplitude (non- uniform array) , and (v) constant grat- 
ing amplitude eio + Aei (uniform array) . 

Figure 4 presents three different types of small and slow 
variations of the grating amplitude on the x-coordinate 
in an initially uniform periodic array. These are linearly 
increasing (figure 4(ii)), decreasing (figure 4(iii)) and si- 
nusoidal (figure 4(iv)) dependences of the grating am- 
plitude with the minimum and maximum values of the 
grating amplitude equal to eio and eio-(-Aei, respectively, 
where Aei <C eio. 

The results of the numerical analysis of the effect of 
small and slow variations of the grating amplitude on the 
scattered wave amplitude inside the array are presented 
in figures 5(i) and 5(n) for Aei = eio/10 — 5 x 10""* 
for two different array widths: (i) L — 12 pm, and (ii) 
L = 20 pm. It can be seen that curves (b) and (c) in 
figure 5(i) coincide with each other. This again reflects 



7 




6 8 
X (|im) 



8 12 
X (|im) 



FIG. 5: The dependences of the relative scattered wave amphtudes on distance from the front array boundary inside the arrays 
(figure 4) with eio = 5 x 10~^, Aei = eio/10, e = 5, do ~ pi/4:, the wavelength in vacuum A = 1pm (the grating has a 
period of 0.58pm and is inclined at an angle of n/8 to the front array boundary); (i) L — 12pm, (ii) L = 20pm. Curves 
(a): uniform arrays with constant grating amplitude eio — figure 4(i). Curves (b): non-uniform arrays with linearly increasing 
ei{x) — figure 4(ii). Curves (c): non-uniform arrays with linearly decreasing ei{x) — figure 4(iii). Curves (d): non- uniform arrays 
with sinusoidal dependence of ei(a;) — figure 4(iv). Curves (e): uniform arrays with constant grating amplitude £io+Aei — figure 
4(v). 



the fact that for narrow arrays, the specific dependence 
of small gradual variations of the grating amplitude does 
not affect the scattered wave amplitudes that arc depen- 
dent only on the average values of the grating amplitude 
in the array. Curve (d) in figure 5(1) differs from curves 
(b) and (c). However, this difference is related not with 
the sinusoidal profile of the dependence of the grating am- 
plitude on the a:;-coordinatc (figure 4(iv)), but with the 
slightly different value for the mean grating amplitude 
in this array. If the mean grating amplitudes were the 
same for all three arrays shown in figures 4(ii)-(iv), then 
curves (b)-(d) in figure 5(i) would have been indistin- 
guishable. Curves (a) and (e) correspond to the uniform 
arrays with the grating amplitudes eio (figure 4(i)) and 
eio + ^ei (figure 4(v)), respectively. 

Figure 5(ii) shows that as the array width increases, 
the shape of the dependence of the grating amplitude on 
the xcoordinate becomes significant. For example, for the 
array width of 20 ym, the scattered wave amplitudes are 
slightly different for arrays with increasing and decreas- 
ing grating amplitude (figures 4(ii) and (iii)) — see curves 
(b) and (c) in figure 5(ii). If the array width is increased 
to 40 or more pm (i.e. L > 37av^), then the scattered 
wave amplitude at any point inside the array is approx- 
imately the same as it would have been in the uniform 
array of the same width and with the grating amplitude 
taken from the point of observation in the non-uniform 
array. That is, the scatteredwave amplitudes are deter- 
mined by local values of the grating amplitude. Note that 
this statement is correct only for slow and small varia- 
tions (imperfections) of the grating amplitude (figure 4). 
It is not valid for arrays where the grating amplitude 
tends to zero somewhere in the array or at its boundaries 



(see figures 2 and 3). If the above statement was valid 
for such arrays, this would have meant that the scattered 
wave amplitude in parts of the array with ei ^ must 
increase to infinity, and due to the diffractional diver- 
gence, it must then be infinite in all other parts of the 
array, which is not possible. 

Increase in Aei results in variations of the scattered 
wave amplitudes that are to a very good approximation 
directly proportional to Aei. Therefore, the graphs in 
figure 5 can be used for evaluation of the scattered wave 
amplitudes for any other value of Aei ei • 



IV. DEAS IN ARRAYS WITH VARYING 
GRATING AMPLITUDE 

As has been seen from the analysis of the previous 
section, magnitude variations of the grating amplitude 
may result in noticeable, but not dramatic, changes in 
the field distribution inside the array. For narrow arrays, 
such changes can even be neglected altogether. However, 
if we consider phase variations in the grating, the situa- 
tion may change very radically, and especially for narrow 
arrays. 

For example, consider three different types of depen- 
dences of the grating amplitude on the .x-coordinate in- 
side a non-uniform array — figures 6(i)-(iii). In figure 
6(i), the magnitude of the grating amplitude is constant 
throughout the array, but the phase of the grating expe- 
riences a stepwise variation </> « 180° in the middle of the 
array, i.e. the sign of the grating amplitude ei (a;) changes 
when crossing the interface between the two sections of 
the array &i x = Li = L/2 — see figure 6(i) and figure 1. 
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In all other parts of the array the phase of the grating is 
not changing (figure 6(i)). Thus we have a non- uniform 
array consisting of two joint uniform arrays with different 
phases of the gratings. In the general case with arbitrary 
stepwise phase shift 4> between the gratings in the joint 
arrays, the dependence of the grating amplitude on the 
a;-coordinate is given by the equations: 



ei{x) 



eio 

eio exp(i^) 




if < a; < L/2, 
if L/2 < X < L, 
otherwise. 



(21) 



If (/) « 180°, then equations (21) give the dependence of 
the grating amplitude presented in figure 6(i). 

The second considered type of non-uniform arrays is 
characterized by a linear dependence of magnitude of the 
grating amplitude on the x-coordinate in such a way, that 
the amplitude of the grating is zero at the front and rear 
array boundaries, as well as at the interface x = Li = 
L/2 (figure 6(ii)). There is also a stepwise variation in 
the phase of the non-uniform grating at x = Li = L/2, 
which is equal to (/) ~ 180°. Thus we have two joint 
non-uniform arrays with linearly varying magnitude of 
the grating amplitude and constant phase of the grating 
in each of the joint arrays (figure 6(ii)). If the stepwise 
variation of the phase between the two joint arrays is 
arbitrary, then the dependence of the grating amplitude 
on the a;-coordinate is given by 



' Aeiox/L 
4eio(L/2 - x)/L 
4eioexp(i(?i)(a;-L/2)/i 
4eioexp(i(/>)(L - x)/L 




if < a; < L/4, 
if i/4 < a; < L/2, 
if i/2 < a; < 3L/4, 
if 3i/4 <x<L, 
otherwise, 

(22) 

where the gradient of changing magnitude of the grating 
amplitude is chosen so that the mean magnitude of the 



^0 



x = L 



(i) constant 




(ii) linear 



(iii) sinusoidal 



FIG. 6: Three different profiles of the grating amplitude with 
a stepwise phase variation <j) ~ 180° at a; = L/2: (i) con- 
stant magnitude of the grating amplitude, (ii) linearly vary- 
ing grating amplitude, and (iii) sinusoidally varying grating 
amplitude. 



grating amplitude in the array described by equations 
(22) is the same as for the array described by equations 
(21). That is. 



/ |ei(a;)|da;= / |ei(a;)|da; = G 

Jo Jl/2 



(23) 



where the constant G is the same for both the depen- 
dences (21) and (22). The reasons for using this condi- 
tion are the same as the reasons for integral (18) being 
the same for the arrays with the grating amplitude given 
by equations (16), (17), and (19). 

The third of the considered arrays with a phase vari- 
ation in the grating is presented in figure 6 (iii). It is 
characterized by sinusoidal variations of magnitude of 
the grating amplitude inside the joint arrays in such a 
way that this magnitude is again zero at the boundaries 
of the periodic array and at the interfaces x = Li = L/2. 
The phase of the grating is again constant inside each 
of the joint arrays, but experiences a stepwise variation 
(j) w 180° ai X — Li = L/2 (where the magnitude of 
the grating amplitude is zero) — figure 6 (iii). If the phase 
variation (j) is arbitrary, then such non-uniform arrays can 
be described by the equations: 

'eio[l - cos(47ra;/L)] if < a; < L/2, 

ei(a;) = I ewe'*{l - cos[27r(2a;/L)]} if L/2<x< L, 
otherwise. 

(24) 

where condition (23) is again satisfied. 




180 185 190 195 
Phase sliift (°) 



200 205 



FIG. 7: The dependences of the relative scattered wave am- 
plitudes at the front boundary on the phase shift (j> for the 
non-uniform arrays described by equations (21), (22), (24) 
with eio = 5 X 10~^, e = 5, 9o = 7r/4, L = 20iim, and the 
wavelength in vacuum A = Ijim. Curve (a): non-uniform ar- 
ray with constant magnitude of the grating amplitude. Curve 
(b): non-uniform array with the linear dependence of ei{x). 
Curve (c) : non- uniform array with the sinusoidal dependence 
of ei{x). Curve (d): uniform array with = 0. 
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FIG. 8: The dependences of the relative scattered wave amphtudes on distance from the front array boundary inside the 
non-uniform arrays described by equations (21), (22), (24) with eio = 5 x 10"'^, e — 5, 9o — tt/4, the wavelength in vacuum 
A = 1pm; (i) and (ii) L = 20pm, (iii) L — 40pm, (iv) L = 80pm. Curves (a): non-uniform arrays with constant magnitude 
of the grating amplitude (see equations (21)); (i) (f> = 190.4° , (iii) = 178.6° , (iv) <j) = 157.6°. Curves (b); non-uniform 
arrays with the linear dependence of ei{x) (see equations (22)); (ii) (f) = 194.7°, (iii) = 182.5°, (iv) (f> = 133.8°. Curves 
(c): non-uniform arrays with the sinusoidal dependence of ei(a;) (see equation (24)); (ii) (p = 195.9°, (iii) (j) = 184.6°, (iv) 
(j) = 129.4°. Curves (d): uniform arrays with (j) — and constant grating amplitude ei. 



The dependences of the ratio of the amplitudes of the 
scattered and incident waves at the front array boundary 
on the phase shift (j) in non-uniform arrays (21), (22) and 
(24) with L — 20 pm and eio = 5 x 10"'^ are presented 
in figure 7 for bulk TE electromagnetic waves that are 
incident onto the arrays at the angle Oq ~ tt/A. This 
figure demonstrates that when the value of the phase 
shift is relatively close to pi, then the amplitude of the 
scattered wave at the front array boundary very strongly 
(resonantly) increases as compared with the conventional 
EAS (cf curves (a)-(c) with curve (d) in figure 7). Thus 
we have two simultaneous resonances in the structure. 
One of these resonances occurs at a resonant frequency 
(wavelength) that is determined by the Bragg condition 
(3). This resonant is obviously common for all types of 
Bragg scattering. In the case of EAS, it results in a strong 
increase of the scattered wave amplitude see curve (d) in 



figure 7. 

The second resonance takes place at an optimal (reso- 
nant) phase shift in a non-uniform array. This resonance 
occurs on the background of an already resonantly large 
scattered wave amplitude that is typical for the conven- 
tional EAS. As a result, the scattered wave ampHtude 
increases many times further as compared with the in- 
cident wave amplitude (see curves (a)-(c) in figure 7). 
This effect was called DEAS }Toj. 

Curve (a) in figure 7 corresponds to the non-uniform 
array with the grating amplitude given by equations (21). 
This is the case of DEAS that was analysed previously 
in 10] for bulk electromagnetic waves. Curves (b) and 
(c) demonstrate that similar strong DEAS also occurs in 
non-uniform arrays with varying magnitude of the grat- 
ing amplitude, e.g. in arrays with the grating amplitude 
given by equations (22) and (24). It can be seen that the 
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maxima of the scattered wave amplitude at the front ar- 
ray boundary are about 15% lower than for curve (a), and 
resonant values of the phase variation are shifted to larger 
values, i.e. to about 195° as compared with 190.4° for 
the array with constant magnitude of the grating — curve 
(a). This demonstrates that slowly increasing grating 
amplitude may result in noticeable changes of resonant 
values of (j) as compared with DEAS in arrays with con- 
stant magnitude of the grating amplitude. On the other 
hand, curves (b) and (c) are very similar, and this indi- 
cates that there is no significant difference in DEAS in 
non-uniform arrays with linearly and sinusoidally vary- 
ing grating amplitude. Thus DEAS proves to be fairly 
insensitive to the particular profile of the slowly varying 
grating amplitude inside the array. 

Physically, the differences between curves (a)-(c) in 
figure 7 are related to the fact that if the grating ampli- 
tude has the profile described by equations (22) and (24) 
(see also figures 6(ii) and 6(iii)), then the amplitude of 
the grating near the middle of the array (where the phase 
variation takes place) is very small, and it is smaller for 
the sinusoidal dependence (24). Thus, the efficiency of 
re-scattering in this area is low. However, this is the re- 
gion near the interface with the phase variation, where 
the diffractional divergence of the scattered wave and its 
re-scattering in the array result in increasing the incident 
wave amplitude, which in turn gives rise to larger scat- 
tered wave amplitudes Therefore, the amplitude of 
the scattered wave must be smaller in arrays with the 
grating amplitude that decreases in the middle of the ar- 
ray (i.e. near the interface with the phase variation in the 
grating) . Thus the resonance with respect to phase varia- 
tion is broader and weaker for curves (b) and (c), than for 
curve (a) in figure 7. Similarly, the resonance for curve 
(c) , corresponding to sinusoidal variations of the grating 
amplitude, is slightly weaker and broader than that for 
the linear dependence (curve (b)). 

The dependences of the scattered wave amplitude on 
the x-coordinate inside the non-uniform arrays (21), (22) 
and (24) of different widths are presented in figure 8. 
Here, the values of the phase shift 4> are chosen so that 
to give maximal scattered wave amplitude at the inter- 
face X = L/2 where the step-like phase variation occurs. 
It can be seen that as the array width decreases from 
40-20 pm, the scattered wave amplitude increases dra- 
matically (more than five times). This feature is typical 
of DEAS in arrays with constant and varying magnitude 
of the grating amplitude (see also O). As can be seen 
from figure 8, L — Lc ~ 40 pm is a critical array width, 
because strong DEAS takes place only ii L < 40iim. 
This critical array width is the same for all the arrays 
(21), (22) and (24), regardless of whether the magnitude 
of the grating amplitude is varying or not. As has al- 
ready been mentioned, the critical array width is deter- 
mined by the distance within which the scattered wave 
may spread along the x-axis before being re-scattered by 
the grating |lCil |. The average grating amplitudes in ar- 
rays (22) and (24) are the same as in array (21) (see 



condition (23)), and therefore this distance is the same 
for the analysed arrays. 

li L > Lc (figure 8(iv)), then the amplitude of the 
scattered wave at the front and rear boundaries of the 
periodic arrays are approximately the same as in the 
case of EAS in the uniform array, because the diffrac- 
tional divergence of the scattered wave from the oppo- 
site side of the array makes only negligible contribution 
to the scattering. However, in the middle of the array 
near the interface with the phase variation, we observe 
noticeable maxima of the scattered wave amplitudes for 
all three considered arrays. The width of these maxima 
is the same and equal (as expected) to the critical array 
width. The maximum corresponding to array (24) is the 
smallest, while the maximum corresponding to the array 
with constant magnitude of the grating amplitude is the 
biggest figure 8(iv). This is again due to the fact that 
the scattering near the interface x = Li in arrays (22) 
and (24) is impeded by the small grating amplitude. 



X = 

I „ 'L/2 



x = 



X = L 



(i) constant 




x = L 



(ii) linear 



(iii) linear 



X = 

I „ '1/2 



x = 



X = L 



X = 
I „ "1/2 



x = 



-e^o-A8i- 



x = L 



(Iv) sinusoidal (v) constant 

FIG. 9: Five different profiles of the grating amplitude inside 
non-uniform periodic arrays. The grating amplitude has small 
gradual magnitude variations, and a stepwise phase variation 
cj) ^ 180° at x = L/2. 

The effect of small and gradual linear and sinusoidal 
variations of magnitude of the grating amplitude (figure 
9) on DEAS is presented in figure 10 for Aei = eio/10 = 
5 X 10^^ and L — 20 pm. Note that unlike EAS in non- 
uniform arrays with small variations of the grating am- 
plitude (figure 5), DEAS appears to be more sensitive 
to the particular shape of the dependence of ei{x) on 
X — cf curves (b) and (c) in figure 10. The sensitivity of 
DEAS to small imperfections of the array is also notice- 
ably stronger than for EAS (the difference between curves 
(e) and (a) in figure 10 reaches up to ~ 50%). This is 
understandable because the resonance during DEAS is 
noticeably sharper (due to the combination of the two si- 
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FIG. 10: The dependences of the relative scattered wave am- 
plitudes at the front array boundary on the phase shift (p 
for the non-uniform arrays presented in figure 9 with eio — 
5 X 10"^, Aei = eio/10, e = 5, 6lo = 7r/4, L = 20 pm, and the 
wavelength in vacuum A = 1 pm. Curves (a)-(e) correspond 
to the non- uniform arrays represented by figures 9(i)-(v), re- 
spectively. 



where d is the thickness of the guiding slab, ^i(a;) is 
the varying grating ampUtude (corrugation amplitude), 
,f{xo) is an arbitrary periodic function with a period 
of 2T:/q, with niax{|/(a;o)|} — 1, and a mean value of 
zero. Dissipation is neglected and all media in contact 
are isotropic. 

For the approximation of slowly varying amplitudes to 
be valid, the corrugation must be small as compared with 
the grating period: 

|6k/(27r)«l. (26) 

This inequality is similar to condition (2) for bulk waves. 

If the corrugation is non-sinusoidal (the function /(xq) 
is not a sine or cosine), it can be expanded into the 
Fourier series 

+ 00 -|-oo 

f{xo)^ ^ fpexp{ipq^x+ipqyy) = ^ fp exp(ipqxo) . 

p—-~00 P— — 00 

In this case, only two complex conjugate terms (similar 
to equations (1)) satisfying the Bragg condition 

ki - ko = -pq (28) 



multaneous resonances) than the resonance during EAS. 

V. GUIDED ELECTROMAGNETIC MODES 

It was mentioned in section 2 that equations (5)-(13) 
and (15) are valid for all types of waves, including guided 
optical modes. This is because the approach based on al- 
lowance for the diffractional divergence of the scattered 
wave, which led to the coupled wave equations (10) and 
(12), is readily applicable for analysis of EAS of guided 
modes in a slab with a periodically corrugated bound- 
ary. However, in this case the coupling coefficients Tq 
and Fi are no longer given by equations (14), but are 
determined in one of the following modern theories of 
mode coupling in corrugated optical waveguides: bound- 
ary perturbation theory Il6|. mode matching theory [l7j , 
local-normal-mode theory Il8|. or direct approximate so- 
lution of the wave equation jT^. These theories are cor- 
rect for any type of polarization (TE and TM) of slab 
modes. Therefore, the results of this paper are also cor- 
rect for any type of guided slab modes. 

The geometry of EAS and DEAS of slab modes is again 
presented by figure 1. In this case the plane of the fig- 
ure is the plane of a slab, and the periodic structure is 
represented by a periodic corrugation of a slab boundary. 
If we consider a non-uniform periodic groove array with 
varying grating amplitude, then the corrugation is given 
by the equations 

U + Ciix)fixo) forO<x<L; ^^5) 
\d for a: < or a; > L, 



must be taken into account. Here, q is the reciprocal 
lattice vector of the grating, which is parallel to the xq- 
axis — see figure 1. 

The vectors ko and ki are the wavevectors of the inci- 
dent and scattered guided modes. Note that in the case 
of guided modes, |ko — |ki| only for TE„-TE„ or TM„- 
TM„ scattering {n is the order of the mode). If EAS is 
related with polarization change (e.g. TE„-TM„ scatter- 
ing), or with change of the mode order (e.g. TE„-TEm 
scattering with n ^ m), then |ko ^ |ki|. That is why 
in equations (8)-(ll) we used ki even though for bulk 
waves fci is equal to ka. 

Therefore, all the speculations, derivations and results 
obtained in sections 3-5 are applicable for EAS and 
DEAS of guided optical modes in periodic groove ar- 
rays with varying corrugation amplitude if ei{x), eio, and 
Aei are replaced by ^ip{x) = £,i{x)fp, ^lop = Cio/p, and 
A^ip = A^i/p, respectively. For example, consider EAS 
or DEAS of an incident zeroth order TE slab mode into 
a scattered zeroth order TE slab mode in the structure: 
vacuum — GaAs slab (with permittivity 12.25) — AlGaAs 
substrate (with permittivity 10.24); a slab thickness of 
d = 0.611m, angle of incidence 6*0 = 7r/4, the wavelength 
in vacuum A = 1.5 pm, and the corrugation is assumed to 
be sinusoidal, i.e. /(xq) = sin(a2:o). If the corrugation 
amplitude in this structure ^1 is given by the equation 
^1 — B.Sei (pm), where ei is the grating amplitude used 
in sections 3-5, then the coupling coefficients Kq and Ki 
are the same as for the bulk TE waves in sections 3-5. 
Therefore, all the graphs from these sections are valid 
for the scattering of the guided slab modes in the above 
structure with the same (as for bulk waves) array widths 
see sections 3-5. 
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VI. CONCLUSIONS 

In this paper, we have analysed theoretically and nu- 
merically the effect of the grating amplitude, slowly vary- 
ing across a non-uniform periodic array, on EAS and 
DEAS of bulk and guided optical waves. The main fea- 
tures of EAS in DEAS in such non-uniform arrays were 
explained by the diffractional divergence of the scattered 
wave. 

In particular, it was shown that the pattern of EAS in 
narrow arrays with gradually varying magnitude of the 
grating amplitude is almost exactly the same as for the 
uniform array with the same width and grating ampli- 
tude equal to the average amplitude of the grating in the 
nonuniform array. On the other hand, in thick arrays, the 
effect of gradual variations in the grating amplitude on 
the scattered field distribution was demonstrated to be 
rather noticeable. At the same time, the main feature of 
EAS — the strong resonant increase in the scattered wave 
amplitude — is typical for all considered non-uniform ar- 
rays. This will allow use of non-uniform periodic arrays 
with gradually increasing grating amplitude for suppres- 
sion of edge effects that are expected to be unusually 
strong during EAS, and for improving the side- lobe struc- 
ture of the scattered signals. 

The tolerance of EAS to small gradual variations (im- 
perfections) in the grating amplitude inside the arrays 
was determined. It was shown that EAS in narrow ar- 
rays is sensitive only to variations of the mean value of 
the grating amplitude, while the particular shape of the 
dependence of the grating amplitude on distance from the 
array boundaries does not matter. In contrast, in wide 
arrays, the scattered wave amplitude has been demon- 
strated to depend mainly on local values of the grating 
amplitude. 

It was also shown that DEAS which was previously 
analysed in arrays with constant magnitude of the grat- 
ing amplitude 10] also occurs in non-uniform arrays with 
varying grating amplitude, though it appears to be more 
sensitive to the particular profile of the grating amplitude 
than EAS. Two strong simultaneous resonances in DEAS 
result in a much greater increase in the scattered wave 
amplitude for the same amplitude of the grating than 
during EAS. This is the reason why edge effects must be 
much more significant for DEAS than for EAS, and the 
use of arrays with gradually varying amplitude is even 
more important in this case. The sensitivity of DEAS to 
small grating imperfections was shown to be noticeably 
stronger than that of EAS, which is related to the much 
stronger resonance taking place during DEAS. 

The approach for the analysis used in this paper, is 
based on allowance for the diffractional divergence of the 
scattered wave. Unlike the dynamic theory of scatter- 
ing IE Hi that is applicable only for the anal- 
ysis of EAS of bulk electromagnetic waves, the approach 
used is directly applicable for the analysis of EAS and 
DEAS of all types of waves, including bulk, guided and 
surface optical and acoustic waves in uniform and non- 



uniform periodic Bragg arrays. Therefore, the obtained 
results are also applicable to EAS and DEAS of Rayleigh 
surface acoustic waves in periodic groove arrays. In this 
case we only need to use in equations (11) and (13) the 
appropriate coupling coefRcients Tq and Fi from the dy- 
namic theory of scattering of Rayleigh waves in periodic 
groove arrays ■ 

The obtained results will be important for development 
of new EAS- and DEAS-based structures and devices in 
optical and acoustic signal-processing, communication, 
instrumentation and sensor design. 
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APPENDIX 

The numerical solution of equations (10) and (12) was 
carried out using the following boundary conditions at 
the array boundaries: 

(dSi/dx),=o = 0, (A.l) 
{dEi/dx),=L = 0, 

where Eqq is the amplitude of the incident wave at a; < 0. 
The scattered wave outside the array is represented (sim- 
ilarly to uniform arrays 0, B 0] by two plane waves 
propagating parallel to the array boundaries one on each 
side of the array (these waves are the solutions to equa- 
tions (15)). Their amplitudes, Ai for a; < 0, and A2 for 
X > L, are determined from two other conditions 

A-i = Ei\x=o+o, A2 = Ei\x=L-o (A. 2) 

that are actually independent of equations (A.l). 

The condition |_Bo|a;=L-o = I^-ooIj representing the 
energy conservation in the steady-state EAS, will be 
satisfied automatically and we do not need to take it 
into account (see also |E S IS Since we neglect 

edge effects at the array boundaries, the conditions im- 
plying the continuity of the derivative dEo/dx across 
these boundaries are not taken into account (this is 
the usual approximation in modern theories of Bragg 

scatteringSlllllillSllllllllilllllllli)- The 

grating amplitude inside the array is assumed to vary ei- 
ther slowly, or very quickly (stepwise variations). These 
variations are taken into account in the function ei{x) in 
the coupling coefficients i^o.i- Therefore, we do not need 
to take into account boundary conditions at interfaces 
with step- wise variations in the grating amplitude inside 
the array (figure 1). 

For convenience, the coupled wave equations (10) and 
(12) are rewritten as a set of three linear first-order dif- 
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ferential equations: 

dEo/dx - iKi{x)Ei{x) = 0, 

dEi/dx - E2ix) = 0, (A.3) 
dE2/dx + Koix)EQ{x) = 0. 

This set of differential equations can be approximated 
by a set of 3 X (iV — 1) finite difference equations (FDEs) 
on a set of N points {xi,X2, ...,xn}- 

Eo{Xn) - En{x,i^i) + — — ^^"-1 [Eo{Xn) + Eo{Xn-l)] 

X [iKi{xn)Ei{x„) +iKi{xn_i)Ei{x„-i)] = 

Ei{Xn) - Ei{Xn-l) + ^'^"-l [Ei{Xn) + Ei{Xn-l)] 

X [E2{Xn) + E2{Xn-l)] = 0, 
Eo{Xn) - Eo{Xn-i) + ^" ^^""^ [Eo{Xn) + EQ{Xn-l)] 
X [-Ko{Xn)E2{Xn) - Ko{Xn-l)E2{Xn-l)] = (A.4) 



where n = 2, 3,...,iV. Since the amphtudcs Eq and Ei 
arc constant outside the array, the points xi and xj^ are 
chosen to correspond to the front and rear boundaries 
of the array, i.e. a;i = and Xn = L. Then boundary 
conditions (A.l) can be written as: 



Eo{xi)=Eoo, E2{xi)=0, E2{xn)=0. (A.5) 



Equations (A.4) and (A.5) form a set of 3A'' 
Unear algebraic equations with 3A^ variables 
Eo{xi),Eo{x2), ...,Eo{xn), Ei{xi),Ei{x2), £^i(xjv), 
E2{xi), -£^2(2^2), E2{xj\[), which was solved numerically 
for different profiles of varying grating amplitude inside 
the array (see sections 3-5). 
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